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ABSTRACT 


An  exponential  relation,  log  N  =  a  -  bM  where  N  is  the  annual 
number  of  earthquakes  of  magnitude  M,  is  generally  accepted  to  exist  between 
the  number  of  earthquakes  occurring  in  a  given  time  and  the  earthquake  magni¬ 
tude.  A  linear  least- squares  solution  of  this  equation  indicates  that  a  small 
earthquake  (2  <  M  <  3)  occurs  every  30  sec.  The  only  difference  between 
these  small  P  waves  and  an  event  is  one  of  magnitude;  therefore,  it  is  reason¬ 
able  to  expect  that  the  P-wave  noise  would  contain  discrete  P-w?.ve  arrivals 
which  travel  across  the  array  as  plane  waves.  The  purpose  of  this  study  is 
to  identify  and  describe  (statistically)  such  discrete  arrivals. 

Two  techniques  have  been  used  to  identify  discrete  arrivals. 

One  method,  Fisher  analysis  of  variance,  provides  an  output  time  function 
which  measures  the  likelihood  of  an  aritval  in  a  specific  segment  of  noise. 

The  other  technique  calculates  a  moving  power  spectra  and  looks  for  spectral 
characteristics  which  are  sufficiently  deviant  from  the  expected  values  of  a 
Gaussian  stationary  process.  Using  an  8-min  noise  sample,  the  methods  were 
able  to  detect  no  discrete  arrivals  propagating  across  the  array  like  plane 
waves.  Within  the  resolving  power  of  these  tests,  no  small  P  wavelets 
appeared  to  be  in  the  mantle  P-wave  nois* 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  report  presents  results  of  an  investigation  of  discrete 
arrivals  in  mantle  P-wave  noise. 

An  exponential  relation  is  generally  accepted  to  exist  between 
the  number  of  earthquakes  occurring  in  a  given  time  and  the  earthquake  mag¬ 
nitude.  Richter  1  gives  this  relation  as  log  N  =  a  -  bM  where  N  is  the  annual 
number  of  earthquakes  of  magnitude  M.  A  linear  least-squares  solution  of 
this  equation,  developed  from  data  of  I960  and  1963, 2  indicates  that  a  small 
earthquake  (2  <  M  <  3)  occurs  every  30  sec.  The  only  difference  between 
these  small  P  waves  and  an  event  is  one  of  magnitude;  therefore,  it  is  reason¬ 
able  to  expect  that  the  P-wave  noise  would  contain  discrete  P-wave  arrivals 
which  travel  across  the  array  as  plane  waves.  The  purpose  of  this  study  is 
to  identify  and  describe  (statistically)  sucn  discrete  arrivals. 

Two  techniques  have  been  used  to  identify  discrete  arrivals. 

One  method,  Fisher  analysis  of  variance,  provides  an  output  time  function 
which  measures  the  likelihood  of  an  arrival  in  a  specific  segment  of  noise. 

The  other  technique  calculates  a  moving  power  spectra  and  looks  for  spectral 
characteristics  which  are  sufficiently  deviant  from  the  expected  values  of  a 
Gaussian  stationary  process.  Theoretical  discussions  of  these  methods  and 
derivations  of  formulas  are  included  in  Section  II. 

Using  an  8-min  noise  sample,  the  methods  were  able  to  detect 
no  discrete  arrivals  propagating  across  the  array  like  plane  waves.  Within 
the  resolving  power  of  the  above  tests  (discussed  in  Section  III),  no  small 
P  wavelets  appeared  to  be  in  the  mantle  P-wave  noise. 
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SECTION  II 

METHODS  OF  DETECTION 


A.  FISHER  ANALYSIS  OF  VARIANCE 

The  analysis -of- variance  technique  is  a  method  of  signal  de¬ 
tection  applicable  only  to  array  data.  This  method  depends  on  the  similarity 
of  signals  across  channels  as  compared  to  the  noise.  The  basic  computation 
provides  an  output  function  of  time,  each  point  of  which  can  be  translated  into 
a  probability  that  a  specific  segment  of  the  original  array  data  contains  a  sig¬ 
nal.  This  technique  has  been  developed  in  previous  literature,^  but  a  brief 
discussion  of  the  theory  is  presented  here  for  the  reader's  convenience. 

The  mathematical  model  for  array  jata  which  has  been  de¬ 
veloped  is 

X. .  =  n  +  t.  +  e  . . 

3  13 


for  i  *  1#  •  •  .  ,  r  and  j  =  1,  .  .  .  ,  s  where  the  parameters  have  the  following 
meaning : 


is  the  observed  value  for  ittl  channel,  jtl1  time 

«  . ,  is  random  noise 
ij 

|i  is  overall  mean  (expected  tr  be  0) 

signal  effect,  alike  across  instruments  in  a  given  subarray 
J  but  different  in  tim>'  and  for  different  subarrays 


The  noise  term  c  is  assumed  to  be  normally  and  independently 
distributed  around  a  mean  of  0  with  an  U'o, known  but  constant  variance 


Thus, 


the  distribution  function  of  any  single  observation  X. .  may  be  written 


11 V 


— ~  exp  -1/2 
av2n 
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and,  from  .e  assumption  of  independence,  the  joint  distribution  of  all  obser¬ 
vations  is 


f  =n  f(x..) 

i.j  y 

With  these  assumptions,  it  is  possible  to  derive  a  criterion 
(Fisher's  F)  to  test  the  hypothesis 

Ho:  t.  =  0 
J 

for  j  =  1,  .  . .  ,  s.  It  can  be  shown  that  F  is  equivalent  to  a  likelihood  ratio 
test  X  of  the  hypothesis  that  there  is  no  similar  signal  effect,  as  compared 
to  the  noise,  across  r  instruments  within  corresponding  noise  segments  of 
time-width  s.  The  likelihood  ratio  test  of  Ho  is  defined  as 


max 

U.O2 

n  f(x  ;  o2,  t  =  0)1 

-i.j  J 

max 

[,n  ft*..;  a2,  tj)j 

[x,aZ.t 

where  max  [L]  is  the  greatest  value  that  the  likelihood  function  L  can  assume. 

It  is  clear  that  \  tends  to  be  less  when  the  ft.'s  are  nonzero  than  when  the  t 's 

J  j 

are  zero,  since  the  likelihood  function  in  the  numerator  cannot  exceed  the 
likelihood  function  in  the  denominator.  Thus,  X  cannot  exceed  1  and  is  always 
positive.  Therefore,  it  is  reasonable  to  reject  the  hypothesis  Ho  when  X  is 
less  than  or  equal  to  some  critical  value  Xc;  i.  e.  ,  the  values  of  X  for  which 
Ho  would  be  rejected  satisfy  0  s  X  *  Xc. 

It  has  been  shown  that  X  and  F  are  related  by  a  monotone  func¬ 
tion  in  such  a  manner  that  large  values  of  F  correspond  to  small  values  of  X; 
thus,  the  critical  rejion  for  the  Fisher  values  is  the  interval  Fc  s  F  s  »  where 
Fc  is  tiie  ci'itical  value  of  F  A  signal  is  considered  present  in  the  noise  if 
F  >  Fc,  and  no  signal  is  present  If  Fs  Fc.  Tables  of  the  distribution  cf  F 

are  readily  available,  and  Fc  can  be  easily  obtained  from  the  desired  frlse- 
alarm  rate.  _ _  _ 
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B.  DYNAMIC  POWER  SPECTRA 


The  dynamic- power-spectra  method  calculates  a  moving  power 
spectra  and  looks  for  anomalously  high  power  in  a  given  frequency  band.  The 
purpose  ia  to  detect  discrete  arrivals  by  looking  for  spectral  characteristics 
which  are  sufficiently  deviant  from  the  expected  values  of  a  Gaussian  stationary 
process 

Consider  a  time  trace  of  length  T,  sampled  at  a  rate  At,  having 

NP  points.  This  trace  is  divided  into  64-  point  overlapping  blocks  so  that  the 

33r^  point  of  one  block  corresponds  to  the  first  point  of  the  subsequent  block. 

The  complex  Fourier  transform  of  each  block  is  taken,  using  the  Cooley- 

Tukey  algorithm.^  Multiplying  the  Fourier  transform  by  its  complex  conjugate 

yields  the  measured  power  spectra  of  the  64-point  data  block.  The  first  32 

positive  frequencies  are  retained,  from  0  to  f  -  A  f ,  in  equal  steps  of  Af 

max 

where 


Af  = 


1 

64At 


and 

"max  2At 

The  average  measured  power  spectrum  of  the  trace  is  found 
by  averaging  the  measured  power  spectra  values  "ta  given  frequency  for  all 
the  nonoverlapping  but  adjacent  data  blocks,  i.  e.  ,  averaging  the  values  at  a 
given  frequency  from  the  first,  third,  fifth,  etc.  ,  data  blocks.  The  measured 
pcv  ^r  spectra  of  the  data  blocks  are  smoothed  over  two  frequencies  and  deci¬ 
mated,  yielding  a  power  spectrum  having  16  frequency  points.  This  method 
is  very  similar  to  one  previously  published  by  Haubrich.  ® 


A  significantly  different  behavior  from  the  average  at  some 
range  in  the  time  trace  may  be  determined  by  noting  if  the  power  spectrum 
of  some  data  block  is  higher  than  the  average  over  some  band  of  frequencies. 


II-3 


•clones  services  division 


This  behavior  is  called  signal.  The  signal  may  be  characterized  by  a  vector 

m  a  p-dimen.ional  space,  each  dimension  corresponding  to  a  frequency  and 
having  components 


Z.  -  N. 
i  1 

a. 

1 

at  the  i,h  frequency  where  2.  ie  the  measured  pow  -  spectrum  (smoothed  or 
unsmoothed),  N.  is  the  average,  and  a.  is  the  standard  deviation.  The  mag- 
nitude  of  the  vector  M  is  then 


where  p  is  the  number  of  dimensions.  Explicit  expressions  for  N.  and  a.  are 
given  by  11 


where  Z. .  is  the  measured  power  spectrum  of  the  jth  data  block  at  the  ith  fre 
quency  and  MB  is  the  number  of  nonoverlapping  data  blocks. 

As  jKnmn  in  Appendix  A,  assuming  that  the  time  trace  is  the 
output  of  a  Gaussian  process,  the  measured  (unsmoothed)  power  spectrum  is 
distributee  as  a  constant  times  a  X2  fund,  -n  having  two  degress  of  freedom, 
and  the  frequency  points  (as  calculated  abo ,.)  are  independent.  The  average 
measured  power  spectrum  is  taken  to  approximate  the  true  power  spectrum. 
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I 


This  approximation  simplifies  the  theoretical  work;  thus,  smoothing  over  two 

2 

frequency  points  corresponds  to  a  X  distribution  with  four  degrees  of  freedom. 


Based  on  experiments  with  traces  of  Gaussian  random  noise 
(autocorrelated),  the  decision  was  made  to  try  two  different  techniques  of 
smoothing  and  analyzing  the  data.  One  consists  of  averaging  the  adjacent  fre- 
quer  cies,  decimating  to  16  frequency  points,  and  looking  at  a  vector  formed 
by  three  adjacent  smoothed  frequencies  (i.  e.  ,  six  original  frequency  incre¬ 
ments).  Only  blocks  having  smoothed  spectral  estimates  which  are  higher 
than  the  average  (N.)  at  all  three  frequencies  are  considered  significant.  Ihe 
magnitude  of  the  resulting  3-dimensional  vector  (M)  is  then  examined  to  give 
greater  discrimination.  The  second  method  would  be  identical  to  the  first  ex¬ 
cept  that  four  adjacent  smoothed  frequencies  (eight  original  frequency  incre¬ 
ments)  are  used. 


Using  these  theoretical  distributions  of  the  spectral  estimates, 
a  fa, Ise -alarm  rate  is  established  for  two  types  of  analysis;  i.  e.  ,  a  probability 
number  that  a  64-point  block  yield  a  frequency-domain  vector  of  given  magnitude 
with  each  component  larger  than  the  average  spectrum  (R)  is  assigned  to  the 
events.  The  details  of  this  derivation,  which  involve  a  simplifying  assumption, 
are  given  in  Appendix  B.  The  results  are  shown  in  Figures  II- 1  and  II-2. 


Figure  II- 1  plots  the  probability  of  occurrence  against  the  mag¬ 
nitude  of  the  3-dimensional  vector.  This  probability  graph  assumes  that  all 


of  the  three  components  exceed  the  estimated  true  spectrum  N.. 


Figure  II-2  plots  the  probability  of  occurrence  vs  magnitude  of 
the  4-dimensional  vector.  This  probability  graph  assumes  that  all  of  the  four 


components  exceed  N,,. 


To  verify  the  reliability  of  the  dynamic-power-spectra  technique 
as  a  method  of  signal  detection,  12  traces  of  autocorrelated  Gaussian  noise  which 
are  4000  points  in  length  and  have  a  sampling  rate  of  0.  1  sec  are  generated. 
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Figure  II-2.  False-Alarm  Rate  for  4-Dimensional  Vector 
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A  sample  average  power  spectrum  of  one  of  these  traces  is 
shown  in  Figure  II- 3 -  An  80-point  seismic  wavelet  having  a  power  spectrum  as 
shown  in  Figure  II-3  and  a  waveform  as  shown  in  Figure  II-4  is  superimposed 
on  the  noise  traces  at  signal-to-noise  ratios  of  0.  7,  1.  1,  and  1.6  for  a  band¬ 
width  of  0.62  to  1.  72  cps  after  smoothing  over  two  points.  This  bandwidth 
consists  of  three  frequency  points.  T1  12  traces  of  4000-point  length  may 
be  considered  as  a  single  48,000-point,  80-min  trace.  Dividing  this  trace 
into  64-point  blocks  implies  that  there  are  750  independent  decisions  as  to 
whether  both  signal  and  noise  are  present  or  whether  there  is  only  noise. 

The  technique  is  then  applied  to  the  data,  and  the  measured 
power  spectra  are  smoothed  over  two  points.  A  signal  is  considered  detected 
in  a  segment  of  noise  if  the  measured  power  spectrum  exceeds  the  average 
over  the  bandwidth  used  and  if  the  magnitude  of  the  resulting  vector  exceeds 
4.  0. 

The  results  are  summarized  below. 

•  There  were  seven  8-sec  signals  having  S/N  =  0.  7; 
the  program  found  two  of  them 

•  There  were  ten  8-sec  signals  having  S/N  =  1.  1; 
the  program  found  five  of  them 

•  There  were  three  8-sec  signals  having  S/N  =  1.6; 
the  program  found  all  of  them 

•  The  program  made  two  false  alarms,  calling 
actual  noise  as  "signal"  only  twice 

The  calculated  false-alarm  probability  is  given  by  2/750  =  0.  0027, 
which  is  in  excellent  agreement  with  the  theoretical  0.0037.  The  program  was 
able  to  detect  the  signals  having  S/N  =1.1  half  of  the  time,  doing  better  than 
the*  0.  27  expected  for  S/N  =  1.  0. 
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Figure  II-3.  Power  Spectra  of  Gaussian  Noise  and 

Signal  Superimposed  on  Gaussian  Noise 
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SECTION  III 

APPLICATION  OF  METHODS  TO  LASA  NOISE 


The  noise  sample  which  was  investigated  for  discrete  arrivals 
was  recorded  at  LASA  on  25  March  1966  from  04:26:  12.  8  to  04:34:  12.  7.  An 
antialiasing  filter  was  applied  to  all  noise  recordings,  and  the  resulting  time 
traces  then  were  decimated.  This  procedure  yielded  data  traces  4800  points 
in  length,  with  a  sample  rate  of  0.  1  sec. 

To  detect  discrete  arrivals  propagating  across  the  large  array 
as  plane  waves,  the  outputs  of  the  A,  B,  C,  and  D  ring  subarrays  were  used. 

B4  subarray  was  dead.  These  12  subarrays  were  chosen  because  they  were 
the  closest  within  the  large  array  (Figure  IH-1).  Any  signal  would  have  had 
similar  characteristics  at  all  of  these  inner  subarrays. 

Reliability  of  the  analysis  of  variance  was  initially  investigated 
using  12  channels  of  generated  random  noise.  The  technique  was  applied  to 
these  traces  with  a  time  window  of  6  sec.  The  maximum  Fisher  value  ob¬ 
tained  was  1.  7,  as  indicated  in  the  output  time  trace  of  Figure  IH-2.  This 
value  was  used  as  the  critical  Fisher  value;  i.  e.  ,  no  signal  was  considered 
present  in  any  segment  of  the  noise  unless  the  Fisher  value  of  that  segment 
was  greater  than  1.  7.  The  actual  critical  value  for  a  failure  rate  of  0.  001, 
as  determined  by  a  table  for  the  parameters  used,  was  1.  67. 

When  this  technique  was  applied  to  12  channels  of  LASA  noise 
data,  it  was  found  that  the  data  should  be  deconvolved  to  obtain  reasonable 
Fisher  values.  By  deconvolving  the  time  traces,  the  seismic  noise  is  more 
nearly  fitted  to  the  model  discussed  in  Section  II;  i.  e.  ,  it  more  nearly  satisfies 
the  assumption  that  the  observations  in  time  are  uncorrelated.  This  operation 
did  not  completely  eliminate  the  correlations  in  time;  therefore,  the  data  did 
not  realize  the  false-alarm  rate  of  0.  001  for  an  F  value  of  1.67. 
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To  use  the  analysis -of- variance  technique  on  ihe  large  array, 
an  output  trace  from  each  of  the  12  inner  subarrays  was  used.  One  sen  set 
of  12  traces  wan  obtained  by  using  a  single  seismometer  output  (trace  1  of 
ring  2,  Figure  III-3)  from  each  of  the  subarrays.  These  traces  were  pre¬ 
whitened  with  a  short  (9 'point  )  deconvolution  filter  designed  on  the  output 
trace  of  subarray  AO.  This  set  of  deconvolved  noise  was  then  beam-steered 
toward  an  area  (NE  China)  which  was  known  to  be  seismically  active  on  25 
March  and  toward  an  area  (S.  Algeria)  which  was  quiet.  LASA  travel-time 
anomalies  were  used  in  calculating  the  amount  of  time  shift. 

The  analysis -of- variance  method  was  applied  to  three  sets  of 
12  traces:  the  data  beam-steered  to  S.  Algeria;  the  data  beam-steered  to 
NE  China;  and  the  nonshifted  deconvolved  data.  A  6- sec  time  window  was 
used.  Figures III-4  through  III- 6  show  the  resulting  Fisher  time  traces.  The 
time  scale  on  all  traces  is  20  sec /in.  ,  and  the  lower  trace  in  each  figure  is 
a  continuation  of  the  upper  trace.  The  dotted  line  represents  the  critical 
Fisher  value  of  1.7. 

The  Fisher  trace  corresponding  to  the  data  beam-steered 
toward  the  quiet  zone  (Figure  III-4)  and  the  Fisher  output  trace  for  the  data 
which  had  not  been  time-aligned  (Figure  III- 5)  had  three  segments  of  noise 
each  with  values  greater  than  the  critical  value.  For  the  output  trace  of  the 
data  beam-steered  toward  the  noisy  area  (Figure  HI-6),  five  segments  had 
values  above  1.  7  Table  IH- 1  lis.s  the  Fisher  values  corresponding  to  each 
of  these  no  as  segments.  Each  letter  in  the  table  represents  a  particular 
segment  of  ti  »  noise,  and  this  letter  is  written  above  that  segment  on  the 
Fisher  output  trace.  This  table  provides  a  means  of  determining  those  noise 
segments  which  have  correspondingly  large  Fisher  values  (>  1.  7)  in  each  of 
the  Fisher  output  traces.  For  example,  noise  segment  J,  the  6-sec  block  of 
noise  255  sec  into  the  trace,  had  a  Fisher  value  greater  than  1.  7  in  each  of 
the  three  traces  (i.  e.  ,  traces  beam-steered  toward  NE  China,  straight  sum¬ 
med,,  and  beam- steered  toward  S.  Algeria). 
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Figure  II1-3.  LASA  Standard  Subarray 


III’ 5 


•elonoe  ••rvlooa  division 


Fisher  Time  Trace  for  Single  Seismometer  Ortput  (No  Time  Shift) 


*  t 

V' 


Seismometer  Outputs  Beam-Steered  toward  NE  China 


••  Segment 


Table  HI- 1 


RESULTS  OF  ANALYSIS-OF- VARIANCE 
AND  DYNAMIC -POWER -SPECTRA  TECHNIQUES 
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These  immediate  results  indicated  that  eight  noise  segments 
contained  possible  discrete  arrivals  since  the  Fisher  values  of  those  seg¬ 
ments  were  greater  than  the  critical  value.  To  determine  whether  these  noise 
segments  contained  discrete  arrivals,  two  other  sets  of  data  from  the.  12  sub¬ 
arrays  were  investigated.  A  discrete  arrival  was  considered  present  in  a 
noise  segment  only  if  that  segment  had  a  Fisher  value  above  1.  7  for  all  three 
sets  beam-steered  toward  the  same  zone. 

A  multichannel  filter  system^*  was  applied  to  the  25  noise  traces 
within  each  subarray.  The  resulting  12  MCF  time  traces  were  prewhitened 
using  a  short  deconvolution  filter  designed  on  the  AO  MCF  output.  As  in  the 
previous  case,  the  traces  were  beam-steered  to  two  locations.  Results  of 
the  analysis -of- variance  technique  are  included  in  Figures  III- 7  through  III-9 
and  in  Table  HI-1. 

In  general,  the  Fisher  traces  obtained  from  the  MCF  outputs 
contained  fewer  noise  segments  with  Fisher  values  above  the  critical  value 
than  did  the  Fisher  traces  obtained  using  single-seismometer  outputs.  Also, 
the  Fisher  values  for  a  particular  MCF  set  were  smaller  than  those  of  the 
corresponding  individual-seismometer  Fisher  trace.  Since  the  MCF  outputs 
were  velocity- limited,  the  large  number  of  high  Fisher  values  from  the  single¬ 
seismometer  output  may  have  been  the  result  of  surface-mode  energy.  The 
noise  data  had  not  been  low-cut-filtered  to  eliminate  the  microseism  peak 
energy. 

The  MCF  Fisher  output  trace  for  the  nonshifted  data  (Figure  IH-7) 
contained  only  one  segment  of  noise  with  a  Fisher  value  of  sufficient  magnitude 
to  indicate  an  arrival.  This  segment,  however,  did  not  correspond  to  any  of 
the  noise  segments  in  the  nonshifted  individual  seismometer  output  trace 
(Figure  IH-5)  which  contained  possible  arrivals.  For  the  data  beam-steered 
toward  S.  Algeria,  there  was  one  agreement  between  peaks  in  the  Fisher 
values  of  the  MCF  output  trace  (Figure  IU-8)  and  the  single-seismometer  trace; 
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this  occurred  for  noise  segment  J,  a  segment  260  sec  from  the  beginning  of 
the  data.  Also,  there  was  only  one  peak  corresponding  in  the  Fisher  traces 
for  the  MCF  traces  and  the  individual  seismometer  traces  beam-steered 

toward  NE  China.  Thus,  two  noise  segments  contained  possible  discrete 
arrivals. 


The  third  set  investigated  by  this  method  was  a  straight-sum 
output  of  the  25  traces  within  each  subarray.  These  12  traces  were  pre¬ 
whitened  and  beam-steered  as  in  the  other  two  cases.  The  analysis-of- 
variance  technique  was  applied  with  a  6- sec  time  window.  Figures  HI-10  and 
H  show  the  Fisher  output  traces. 


There  were  seven  segments  of  noise  with  Fisher  values  above 
1.  7  beam-steered  to  NE  China.  For  the  straight  sums  beam-steered  toward 
S.  Algeria  (the  quiet  zone),  only  one  segment  of  noise  had  a  value  (1.  77) 
greater  than  the  critical  value.  There  was  no  agreement  between  this  segment 
(of  the  S.  Algerian  data)  and  noise  segment  J  which  had  a  large  Fisher  vzlue 
in  the  corresponding  individual- seismometer  output  and  the  MCF  output. 
However,  when  the  Fisher  traces  were  investigated  for  discrete  arrivals  by 
comparing  the  three  traces  which  had  been  beam-steered  toward  NE  China, 
one  segment  was  found  to  consistently  have  large  F  values.  For  these  data 
(Figures  IH-6,  III-9.  and  HI-H),  one  segment  275  sec  into  the  record  had  a 
peak  F  value  when  the  values  were  calculated  using  the  individual -seismometer, 
straight-sum,  and  MCF  outputs.  Thus,  the  analysis -of- variance  technique 
which  was  applied  to  a  portion  of  the  large  array  detected  a  possible  discrete 

mantle  P-wave  arrival  in  the  noise  when  the  data  were  beam- steered  toward 
NE  China. 


The  dynamic-power-spectra  technique  was  also  used  to  try  to 
detect  discrete  arrivals  traveling  across  the  large  array  as  plane  waves. 
This  method  was  applied  to  five  single  traces  obtained  by  summing  the  12 
traces  within  each  of  five  sets  used  in  the  analysis -of- variance  study. 
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Fisher  Time  Trace  of  MCF  Outputs  Beam-Steered  toward  S.  Algeria 


of  MCF  Outputs  Beam-Steered  toward  NE  Chin, 


- 0.313  -  I. DM  CPS  BANDWIDTH 

- 0.625  -  1. 06  CPS  BANDWIDTH 

Fisher  Time  Tree,  of  Straight-Sum  Outputs  Beam-Steered  to  S.  Algeria 


The  MCF  outputs  were  not  cammed.  Specifically,  three  traces 
for  the  individual-seismometer  outputs  (NE  China,  S.  Algeria,  no  shift)  and 
two  traces  for  the  straight-sum  outputs  (NE  China,  S.  Algeria)  were  used. 

The  measured  power  spectra  were  smoothed  over  two  frequency  points  and 
investigated  in  the  frequency  bandwidths  of  0.313  to  1.094  cps  and  0.625  to 
1.406  cps.  These  frequency  ranges  correspond  to  a  3-dimensional  vector 
having  a  false-alarm  probability  curve  as  shown  in  Figure  II- 1.  Results  of 
this  method  have  been  superimposed  on  the  Fisher  output  trace  of  Figures  III-4 
through  III- 6  and  Figures  III- 10  and  III- 11.  The  bars  below  the  traces  indicate 
the  64-point  noise  segments  with  measured  power  spectrum  exceeding  the 
average  powei  spectrum  over  the  bandwidths  used.  The  magnitude  of  the 
vector  is  shown  above  each  bar. 

According  to  iigure  II- 1,  the  probability  of  occurrence  of  a 
3-dimensional  vector  with  a  magnitude  greater  than  4.0  is  less  than  0.005. 
Initial  results  of  this  method  when  applied  to  L.ASA  noise  over  the  large  array 
indicate  that  a  reliable  analysis  can  be  made  if  vectors  with  magnitudes  of 
2.5  are  considered.  The  corresponding  probability  of  occurrence  for  3- 
dimensional  vectors  with  magnitudes  above  this  is  less  than  0.019.  There¬ 
fore,  analysis  of  the  dynamic-power-spectra  method  is  based  mainly  on  vector 
magnitudes  greater  than  2.  £ 

The  dynamic  -power  results  for  the  individual-seismometer 
outputs  beam-steered  to  S.  Algeria  (Figure  III-4'  contained  four  segments  of 
noise  with  vector  magnitudes  exceeding  2.  5  for  tne  two  bandwidths  considered. 
For  the  bandwidth  of  0.313  to  1.094  cps,  one  segment  of  noise  with  magnitude 
3.  55  corresponded  to  a  peak  in  the  Fisher  output  trace.  The  other  three  seg¬ 
ments,  however,  did  not  correspond  to  any  peaks  with  Fisher  values  above 
1.  7.  In  particular,  the  measured  spectra  of  noise  segment  J,  at  259  sec  with 
a  maximum  Fisher  value  of  2.  2,  did  not  exceed  the  average  spectra  for  the 
frequencies  within  the  two  bandwidths. 
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For  the  individual-seismometer  outputs  which  were  not  beam- 
steered,  this  method  resulted  in  six  magnitudes  above  2,5.  One  segment 
with  a  Fisher  value  of  2.0  had  a  corresponding  vector  of  magnitudes  3.  7  and 
2.6  for  the  bandwidths  of  0.313  to  1.094  cps  and  0.625  to  1.406  cps,  respec¬ 
tively.  There  was  no  agreement  between  the  other  peaks  and  vectors.  These 
results  are  illustrated  in  Tabic  III- 1 .  The  third  trace  (single-seismometer 
outputs  beam-steered  toward  NE  China)  had  11  magnitudes  greater  th?n  2.5; 
seven  of  these  vectors  corresponded  to  Fisher  values  above  the  critical  value. 
The  segment  of  noise  having  the  largest  Fisher  value  also  had  the  largest 
vector  magnitudes  in  both  bandwidths.  Noise  segment  K  had  four  correspond¬ 
ing  vectors  with  magnitudes  above  2.4.  In  all  cases,  most  vectors  of  large 
magnitudes  (>  2.0)  corresponded  to  peaks  in  the  Fisher  time  traces,  although 
the  values  of  these  peaks  were  below  1.7. 

The  dynamic-power-spectra  technique  was  applied  also  to  the 
sum  of  the  straight  sums  of  the  individual  seismometers  beam-steered  to 
NE  China  and  to  the  sum  of  the  titra  ;ht  sums  beam- steered  to  S.  Algeria. 

For  the  NE  China  data,  three  segments  of  noise  for  the  bandwidth  of  0.313  t 
1,094  cps  and  three  segments  for  the  bandwidth  of  0.625  to  1.406  cps  had 
magnitudes  above  2.5.  One  of  these  segments,  which  had  a  vector  magnitude 
above  2.  5  for  both  bandwidths,  corresponded  to  noise  segment  K  having  a 
Fisher  value  of  1.  76.  For  the  data  beam-steered  toward  S.  Algeria,  there 
were  three  vectors  of  magnitudes  above  2.5  in  the  bandwidth  of  0.313  to  1.094 
cps  and  two  in  the  bandwidth  of  0.  625  to  1.406  cps.  Of  these,  only  one  (with 
a  magnitude  of  3, 5)  corresponded  to  a  Fisher  value  above  the  maximum  value 
of  the  random  noise  data. 

Cuinparing  the  results  of  the  dynamic-power-spectra  method 
for  the  three  sets  of  output  traces,  only  one  consistent  agreement  between  the 
location  of  the  arrivals  was  noted.  In  other  words,  for  one  segment  of  noise, 


ni-18 


•oloneo  ssnrlon  division 


the  power  spectra  exceeded  the  average  power  a,  the  three  frequencies  within 

the  two  bandwidths  and  the  vector  magnitudes  were  greater  than  2.  5  for  both 

the  individual-seismometer  outputs  and  the  straiakt  0 

p  na  tne  straight- sum  outputs.  This  par¬ 
ticular  segment  occurred  in  the  NE  China  data. 

It  was  expected  that,  if  a  discrete  arrival  were  detected  whin 
the  same  segment  of  noise  by  both  methods,  an  arrival  wouid  be  present  in 
the  notse.  This  occurred  only  for  noise  segment  K  for  the  data  beam-steered 

FlSher  Values  were  abore  *•  7  for  this  segment  in  the  MCF, 
straight- sum,  and  single-seismometer  outputs,  and  the  magnitudes  of  the 
vectors  corresponding  to  this  segment  were  greater  than  2.  5.  Thus,  the 
large-array  study  indicated  a  discrete  arrivai  from  NE  China,  behaving  as 
would  be  expected  of  a  small  signal  from  teleseismic  distances. 

In  investigating  the  large-array  outputs,  the  wave  velocity  was 
a  critical  factor  since  the  maximum  separation  between  the  subarray,  was 
54  km.  Thus,  only  a  very  limited  region  of  wavenumber  space  could  be  in¬ 
vestigated  for  discrete  P-wave  arrivals  by  any  one  se,  of  time  shifts.  I„ 
order  to  have  general  simultaneous  coverage  of  the  high-velocity  field,  the 
two  detection  techniques  were  applied  to  output  traces  of  individual  subarrays. 

he  seismometers  on  the  four  inner  rings  (Figure  III-3)  were  used  to  make 
t  e  method,  less  velocity-sensitive.  Only  the  outputs  of  the  four  subarray, 

.  Bl.  S2.  and  B3  were  studied.  To  remove  the  effects  of  microseisms  a, 
ower  frequencies,  the  noise  recordings  were  bandpass-filtered  (0.8-  to 
5.  0-cps  bandpass).  These  filtered  outputs  were  then  prewhitened  using  a 
9-point  deconvolution  filter  designed  on  trace  21  of  each  subarray. 

To  estimate  the  power  of  both  methods  for  the  detection  of 
discrete  arrivals  within  a  subarray,  a  5-.ec  signal  which  had  been  bandpass- 
filtered  and  deconvolved  was  superimposed  on  the  13  filtered  inner  traces  of 
subarray  Bl.  The  same  signal  waveform  was  added  to  each  trace  (filtered 
signal  shown  in  Figure  111-12).  The  signal  was  buried  within  the  8-min  noise 
sample  at  six  different  positions  corresponding  to  six  velocities  of  propagation. 
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The  signal  was  superimposed  on  the  noise  with  signal-to-noise  ratios  of 
approximately  1  and  1/2  To  do  this,  the  signal  was  scaled  by  a  factor  de¬ 
termined  by  dividing  the  maximum  value  of  the  noise  segment  by  the  maximum 
value  of  the  signal.  Figure  IH- 13  shows  a  segment  of  the  time  traces  con¬ 
taining  the  signal  with  an  S/N  =  1/2  and  a  velocity  of  14  km/sec. 

The  Fisher  analysis -of-variance  method,  using  a  time  window 
of  3  sec,  was  applied  to  the  13  traces  containing  the  embedded  signals.  The 
output  trace  is  shown  in  Figure  III- 14.  This  method  was  able  to  detect  all 
signals  with  o/N  -  1/2  o-  larger  and  a  velocity  greater  than  10  km/sec.  The 
dotted  line  across  the  trace  in  Figure  IH-14  corresponds  to  the  maximum  F 
value  obtained  when  the  analysis  of  variance  was  applied  to  the  13  traces  be¬ 
fore  the  signals  were,  superimposed.  Fisher  values  for  the  segments  con¬ 
taining  a  signal  ranged  from  4.  0  to  12.  17. 

The  dynamic-power-epeutra  technique  also  was  applied  to  one 
of  the  traces  containing  the  embedded  signal.  Three  bandwidths  -  0.625  to 
1.406  cps,  1.25  to  2.031  cps,  ar.4  0.  938  to  1.  719  cps- were  used;  each  con¬ 
tained  three  frequency  points.  The  outputs  of  this  method  are  superimposed 
on  the  Fisher  output  time  trace  in  Figure  III- 14.  The  number  of  successful 
detections  varied  among  the  three  bandwidths,  but  the  1.25-  to  2.031-cps 
range  detected  move  of  the  signals.  There  were  only  two  false  alarms:  one 
in  the  0.  625-  to  1.  406-cps  bandwidth  with  a  magnitude  of  2.  7  and  the  other  in 
the  1.  25-  to  2.  031-cps  bandwidth  with  a  magnitude  of  1.  13.  (The  latter  false 
alarm  is  not  shown  in  the  figure  since  it  occurred  on  the  last  part  of  the  trace 
away  from  a  signal. )  These  two  magnitudes  corresponded  to  false-alarm 
probabilities  of  0.  018  and  0.  044.  It  should  be  noted  that  the  results  of  the 
dynamic-power-spectra  technique  were  not  affected  by  the  wave  velocity  of  the 
signals  in  this  experimental  case  since  only  one  of  the  individual  traces  was 
used.  Therefore,  the  results  of  this  method  on  the  experimental  signal  were 
studied  for  the  number  of  signals  detected  and  the  signal-to-noise  ratios. 
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Figure  *11-13.  Signal  Superimposed  on  Noise  of  Subarray  B1 
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Trace  of  Noise  Containing  Embedded  Signals 


This  test  indicated  that  both  the  analysis-of- variance  method 

and  the  dynamic-power-*pectra  technique  should  detect  discrete  P-wave  ar- 

rivals  which  have  a  A  >  30°  and  a  signal-to-noise  ratio  of  1/2  or  larger  and 
are  highly  coherent. 

The  analysis  -  of- variance  technique  was  applied  to  the  13  fil¬ 
tered  inner  traces  of  the  four  subarrays  for  the  .an,,  noise  used  in  the  large- 

array  analysis.  A  3-sec  time  window  was  used.  Resulting  Fisher  time 

traces  are  shown  Figures  111-15  through  HI- 18.  For  subarray  AO.  eight 
segments  of  noise  had  Fisher  values  above  3.0.  The  maximum  value  was 
3.  83.  which  was  below  the  minimum  value  (4.  0)  obtained  in  the  test  case  for 
a  signal  traveling  at  1 1  km/sec  with  S/N  =  1/2.  ln  subarray  B1  shown  i„ 
Figure  HI-18,  four  segments  of  noise  had  Fisher  values  between  4.  0  and  4.  74 
These  values,  compared  to  the  results  of  the  embedded  signal  output  trace, 
indicated  four  possible  discrete  arrivals.  There  were  also  four  Fisher  value, 
urger  than  4. 0  the  output  trace  of  subarray  B2  (Figure  IH-17),  although 
the  largest  value  was  only  4.  15.  Two  possible  discrete  arrival,  were  de¬ 
tected  in  subarray  B3.  corresponding  to  Fisher  value,  of  4.4  and  9.09  (Fig¬ 
ure  IH- 18).  The  largest  value  should  be  highly  significant  and  should,  in 
feet,  be  visible  in  the  time  traces  of  B3.  Figure  111-19  contain,  the  segment 
of  the  13  filtered  traces  corresponding  to  this  large  Fisher  value.  The  signal 
was  clearly  present.  A  small  spike  appears  to  have  caused  this  Urge  value. 

If  a  discrete  arrival  were  present  in  the  noise,  the  Fisher 
time  traces  should  have  a  large  value  for  this  segment  in  all  four  subarray 
outputs.  Therefore,  to  determine  if  the  peaks  in  the  Fisher  output  traces 
corresponded  to  the  same  segments  of  noise  within  each  of  the  subarrays  a 
straight  sum  of  the  four  Fisher  traces  was  calculated.  This  trace  is  shown 
in  Figure  UI-20.  The  sum  indicated  only  two  possibilities  of  detection. 
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Figure  III- 19.  Noise  Segment  Containing  Signal  Detected  by  Analysis- 
_ _ of-Varianee  Technique  Applied  to  Subarray  B3  _ 
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The  maximum  Fisher  values  of  the  first  segment  were  2.  68, 
4.74,  3.42,  and  4.37  for  subarrays  AO,  Bl,  B2,  and  B3,  respectively.  An 
inspection  of  the  13  filtered  traces  for  each  subarray  in  this  area  of  the  data 
revealed  no  visible  signal.  The  second  segment  of  noise  had  Fisher  values 
of  2.3,  4.62,  2.3,  and  9.09  for  the  four  subarrays.  A  signal  was  visible  in 
the  filtered  time  traces  of  Bl  and  B3. 

To  determine  the  effect  of  wave  velocity  on  the  Fisher  outputs, 
the  four  Fisher  traces  were  time-shifted  to  correct  for  the  moveout  of  a  sig¬ 
nal  propagating  with  a  velocity  of  10  km/sec  across  the  four  subarrays  from 
the  north.  Then,  the  resulting  traces  were  summed.  This  procedure  was 
repeated  for  signals  traveling  with  the  same  velocity  from  the  NW,  SW,  S, 

SE,  and  WE.  Time  corrections  were  also  made  for  a  signal  propagating  at 
16  km/sec  from  the  six  directions.  The  resulting  12  traces  had  the  same 
two  predominant  peaks  as  the  sum  of  the  Fisher  output  traces  which  had  not 
been  shifted  (Figure  III- 20). 

The  dynamic-power-spectra  technique  was  applied  to  the  mul¬ 
tichannel-filtered  outputs  of  AO,  Bl,  B2,  and  B3  in  the  three  frequency  band- 
widths  used  in  the  embedded  signal  test  case.  The  64-point  segments  of  noise 
which  this  method  found  to  contain  signals  are  indicated  on  the  Fisher  outputs 
m  Figure  FI-15  through  III-18.  In  each  of  these  cases,  most  noise  segments 
with  vector  magnitudes  greater  than  2.5  corresponded  to  peaks  in  the  Fisher 
time  traces.  However,  the  noise  segments  in  the  Fisher  outputs  having  the 
largest  Fisher  values  did  not  have  corresponding  vectors  of  large  magnitude; 
also,  the  dynamic-power-spectra  method  did  not  detect  arrivals  occurring  at 
the  same  time  within  the  four  arrays. 

In  general,  for  both  methods,  there  was  a  lack  of  correlation 
among  the  four  subarrays  of  the  noise  segments  found  to  contain  possible  dis¬ 
crete  arrivals.  In  particular,  neither  technique  indicated  that  an  arrival  was 
present  in  the  segment  of  noise  which  was  found  to  contain  an  arrival  by  the 
large-array  analysis. 
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It  appears  unlikely  that  any  highly  coherent  small  P-wave 
arrivals  were  present  in  the  noise  sample  studied.  A  single  discrete  arrival 
was  indicated  from  the  NE  China  area,  but  the  "event"  did  not  show  up  in  the 
subarray  processing.  This  could  have  been  masked  by  poorer  S/N  ratios  for 
the  subarray  processing.  Within  our  ability  to  discriminate,  therefore,  it 
appears  that  the  P-wave  noise  does  not  contain  a  significant  component  com¬ 
posed  of  discrete,  coherent  wavelets.  Modeling  of  the  P-wave  noise  as  a 
stationary  Gaussian  process  appears  to  be  most  reasonable. 
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APPENDIX  A 


X  DISTRIBUTION 


Consider  a  time  trace  f{t)  which  has  a  normal  distribution 


[N(0,  o  )] .  The  Fourier  transform  yields 


i.  =  *  (t  )  cos  (2rrkn/N)  k  =  0,  1,  ....  N  -  1 

k  n=0  n 


b  =  53  I  (t  )  sin  (2rrkn/N) 
k  n=C  n 


If  the  f(t  )  are  independent,  a  is  distributed  as  N  (0,  A)  where 

^  K 

2  N“ 1  7  z\  N_1  1  2 

A  =  a  53  coa  {2nkn/N)  =  a  N/2  +  53  ccs  (4rrkn/N)  =  o 


for  k  4  0  or  N/2  since 


4rrkni 


N-l  N- 1  4rrki  _  . 

J3  cos  4nkn/N  =  real  53  e  =  real  --  ; ' 

n*0  n=C  5^- 

e  N  -1 

for  k  ^  0  or  N/2. 


Similarly,  b^  is  distributed  as  N(0,  B)  where 


=  0  +  Oi 


2  N~ 1  2  2  N  N_1 

B  =  a  53  8*n  (2nkn/N)  =  a  —  -  53  sin  (4rrkn/N) 

n=0  n=C 


for  k  /  0  or  N/2. 
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Thus,  A  =  B  and  a  and  b  are  distributed  as  N(0,o'  )  and 

JK  JK 

a  /o'  and  b  /o'  are  distributed  as  N(0,  1).  Squares  of  standard  normal  vari- 
k  k  2  ,2 

ables  are  distributed  as  \  with  one  degree  of  freedom.  Thus,  a^  o  + 

b^/n'^  is  distributed  as  yp  {two  degrees  of  freedom),  or  a2  +  b?  is  distributed 
k  y  K 

as  o,2x|* 

The  independence  of  adjacent  points  in  frequency  is  shown  by 


for  j  4  k.  Using  the  independence  of  f(t  )  f(t  )  with  m  4  n  and  since  the  cos 

n  m 

7 

and  sin  functions  are  orthogonal  under  this  type  of  summation,  one  has 

N-l 

cov  (a.a,  )  =  o  5^  cos  (^Ttjn/N)  sin  (2TTkn/N)  =  0 
J  k  n=0 

Q 

for  j  4  k.  These  results  are  true  asymptotically  as  N  becomes  large. 
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DERIVATION  OF  FALSE-ALARM  RATE 


The  functional  form  of 


2 

*4 


is 


x 

4 


-x/2 

e 


for  x  a  0  and  0  for  x  <  0 


This  has  a  mean  of  4  and  a  variance  of  8. 


The  frequencies  calculated  are  independent  (Appendix  A).  Thus, 

an  n-dimensional  frequency  vector  of  the  type  used  in  this  analysis  would  have 

2 

components  which  are  independent  and  distributed  ay  .. 

A4 

First,  consider  a  2-dimensional  case  which  is  easy  to  visualize; 
then,  generalize  from  this  case. 


Given  a  vector  y  =  (x^,  x^)  where  x^  and  x^  are  distributed 
ax4  and  independent,  we  wish  to  know  the  probability  that  x  >  N.,  x~  >  N~, 

X  1  b  L 

and  magnitude 


where 


MB  (X.  -  N.) 

E  — - 


MB  -  1 
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Assuming  the  estimates  of  the  mean  and  the  standard  deviation 

to  be  correct,  this  is  equivalent  to  wishing  to  know  the  probability,  if  given 

2 

a  vector  y  =  (x^,  x^)  where  x^  and  x^  are  distributed  and  independent,  that 
Xj  >4,  >  an<* 


This  is  the  integral  of  the  joint  distribution  of  x.x_  over  the 

X  b> 

shaded  area  shown  below. 


Since  the  components  are  independent,  the  joint  distribution 
is  the  product  of  the  individual  distributions.  One  then  has  an  integral  of  the 
form  J j  ~  e  ^  e  dx  dy  to  evaluate. 

A1 

When  the  preceding  is  extended  into  three  and  four  dimensions, 
the  integral  becomes  quite  complex.  After  a  short  attempt  to  get  a  closed- 
form  expression  for  these  integrals,  the  attempt  was  abandoned  and  a  sim¬ 
plifying  assumption  made. 
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The  circular  sector  (in  two  dimensions)  is  replaced  by  a  square 
of  equal  area.  To  give  an  upper  and  lower  bound,  the  circular  sector  is  re¬ 
placed  by  an  inscribed  square  and  a  circumscribed  square.  (See  figure  below.  ) 


Extending  this  reasoning  to  higher  dimensions,  the  probability 
integral  for  the  3 -component  (dimensional)  case  becomes 


and  for  the  4- component  (dimensional)  case. 


where  a  is  the  length  of  a  side  of  the  cube  or  hypercube. 

These  integrals  were  evaluated  as  a  function  of  R.  The  re¬ 
sults  are  shown  in  Figures  II- 1  and  II- 2. 

Experience  with  generated  Gaussian  noise  suggests  that  the 
equal  "volume”  curves  give  good  agreement  with  experimental  results. 
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